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Abstract 

' 

| New elements of the dual cone of the set of fermion iV-representable 

2-density operators are proposed. So far, the explicit form of the cor- 
responding necessary conditions for TV-representability is obtained for 
N = 3. In this case the new condition is stronger than the known B- 
i-pH | and C-conditions for 3-representability. The results provide evidence 

that in the spectral decomposition of the iV-representable 2-density 
operator there exists an intrinsic relation between the eigenvalue and 
the corresponding eigenfunction. 
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l. Introduction 



The iV-representability problem [1, 5, 13, 18, 19, 23, 25] appears in non- 
relativistic quantum mechanics of iV-fermion systems. The Hamiltonian 
H N for systems of iV-fermions contain only operators with 1- and 2-body 
interactions: 

N 

i=l l<i<j<N l<i<j<N 
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The ground state energy of the system can be determined variationally by 
minimizing the iV-particle functional 



E = inf Tr(H N D N ) 

over the set V N consisting of all fermion iV-particle density operators. In- 
stead, because of the appearance of at most 2-body interactions between the 
particles, the ground state energy of the system of 7V-fermions could be, in 
principle, determined variationally by minimizing the 2-particle functional 

E= inf (^)Tr(h 2 D 2 ) 



D 2 eV 2 



N 



over the set Vjj consisting of all fermion 2-particle reduced density operators, 
i.e. such 2-particle density operators D 2 which possess an iV-particle fermion 
preimage D N : 

D 2 = L 2 N D N = Tr, X I)\ 

It is known [1, 19] that the set Vjj is a proper convex subset of V 2 , the set of 
all fermion 2-density operators. However, the complete characterization of 
V\j as a convex proper subset of V 2 is not yet known. It has been shown [19] 
that a knowledge of all exposed points (which are extreme) of V% is sufficient 
to characterize the closure of Vfj\ only some of them are known. The dual 
characterization of Vfj involves a determination of the dual (polar) cone 
consisting of 2-particle self-adjoint operators b 2 for which Tr(b 2 D 2 ) > 0, 
\/D 2 e Vf{, that is equivalent to the positive-semidefiniteness of the iV- 
particle operator T^b 2 = b 2 A I^( N ~ 2 ) = A N b 2 ® I®( N ~ 2 )A N > 0. The 
dual cone Vfj is a convex one. Any element of the dual cone Vfj provides 
an iV-representability condition. Those coming from the extreme elements 
of Vfj are the strongest ones. They give the hyperplane characterization 
of and thus the solution of the iV-representability problem. Several 
necessary conditions for A^-representability have been derived and some of 
their structural features and mutual interrelations are established [2, 3, 6, 
7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 19, 20, 21, 22, 24, 26, 27]. 

In this paper we find new elements of the dual cone Pfj, and thus define 
new necessary conditions for iV-representability of a trial 2-particle density 
matrix D 2 . In the case N=3, in which the new elements of the dual cone 
are given explicitly, they lead to necessary conditions for 3-representability 
of the trial D 2 which are stronger in comparison with the known B and 
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C conditions [3, 4, 21]. The results show that in the spectral decomposi- 
tion of a 2-particle fermion density operator D 2 there exists an intimate 
relation between the eigenvalue and the corresponding eigenfunction which 
has to be satisfied in order that D 2 can be 3-representable ( in general N- 
representable). More precisely, the condition obtained for 3-representability 
says that the upper bound of the eigenvalue of D 2 is a functional of the 
corresponding eigenfunction. It is worthwhile to remember that for N- 
representability of a 1-particle density matrix D 1 the upper bound on the 
eigenvalues does not depend on the eigenfunctions. This explains why the 
iV-representability problem for D 2 is so much harder than for D 1 . 



2. The dual P-condition. 

In this paper the underlying 1-particle Hilbert space H 1 is finite dimen- 
sional dimTL 1 = n. TL A2 is a 2-particle antisymmetric Hilbert space, the 
Grassmann product TL 1 A TL 1 = A 2 ^ 1 <8> TL 1 . P 2 = g 2 ® g 2 is the projection 
operator onto a 2-particle antisymmetric function g 2 6 TL A2 . The opera- 
tor ( N 2 )P 2 A J^ 2 ) = ( N 2 )A N P 2 ® I*(»-*)A» = A N Ei<^ Pftj) ® 
I N ~ 2 (1, . . . ,i— • • • ,j — 1, J + l, • • • ,N)A N is the simplest (elementary) 
antisymmetric operator with "2-body interactions" acting on 7i AN . Because 
(^)-^g A I A ( N ~ 2 ) is positive-semidefinite, P 2 belongs to the dual cone V^, 
and therefore Tr(D 2 P 2 ) > is a necessary condition for iV-representability 
of D 2 (the P-condition). Besides, there exists another element of V\ that is 
generated by the 2-particle antisymmetric function g 2 £ 7i A2 . Let A max (g 2 ) 
denote the maximal eigenvalue of {^)Pg Al A ( N ~ 2 \ which in general depends 
on g 2 , then the operator A max (g 2 )I AN - ( N 2 )P 2 A I A ( N ~ 2 ) > (is positive- 
semidefinite), and therefore the 2-particle operator A "7aj [ 9 - I A2 — Pq belongs 
to the dual cone Vfj . Thus, we have obtained 
Theorem 2.1: For any g 2 G H A2 the operator 



Amaxjg 2 ) T /\2 p2 c -52 / 911 
JNj 1 ~ ^9 G ' W> 



cmo? gives i/ie following necessary condition for N- representability of a 2- 

nax (g 



partade fermion density matrix D 2 : Tr[(^f^-I A2 - Pg)D 2 ] > 0, i.e. 



Tr(D 2 P 2 ) < Ama ^ (g ; , V 9 2 eW A2 (2.2) 
(.2 ) 
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Here, A max (g 2 ) is the maximal eigenvalue of the operator {^)Pg f\ I A ( N 2 \ 

We propose to call the new condition "the dual P-condition" . 
In particular, from the above theorem follows 

<") 

Theorem 2.2: If D 2 = Yl}=i ^i^g t ^ s the spectral decomposition of D 2 , 
then it is a necessary condition for N -representability that the eigenvalues 
Aj must satisfy the inequalities \ < Am ^fr ^ . Here, A max (gf) is the maximal 

I 2 J 

eigenvalue of the operator (^P^ A I A ( N ~ 2 ) ; where gf is the eigenfunction 
corresponding to the eigenvalue Aj. 

This theorem shows that the bound on the eigenvalue is a functional 
of the corresponding eigenfunction. For the iV-fermion 1-body elementary 
operator NPg A i 7 ^ -1 ), g 1 e H 1 , the maximal eigenvalue A max (g 1 ) = 1 
for any g 1 , because NPg A / A ( Ar_1 ) is a projection operator. Hence, the 
equivalent of Theorem 2.1 says that Tr(D 1 Pg) < for arbitrary g 1 G 
H A1 , which in particular means that the bound on an eigenvalue of an 
iV-representable 1-density operator does not depend on the corresponding 
eigenfunction, contrary to the 2-density matrix case. 

As seen from Theorem 2.1, the new necessary condition for iV-representability 
of a 2-particle density matrix D 2 requires knowledge of the maximal eigen- 
value A max (g 2 ) of the N- particle operator (^)P 2 A I A ( N ~ 2 ). It is rather 
hopeless to find A max (g 2 ) for arbitrary N and any g 2 G 7i A2 . But even the 
solution for some "simple" g 2 would contribute to knowledge of the struc- 
ture of both the dual cone V\ and the set V\ of fermion A^-representable 
2-density operators. So far, we have results for N = 3, in which case it 
was possible to find the spectral decomposition of the elementary operator 
3P 2 A I 1 for arbitrary g 2 £ Tt A2 . The details and proofs concerning this 
spectral decomposition and the reduced 2- and 1-particle density operators 
corresponding to the eigenstates of the operator 3P 2 A I 1 will be published 
in a separate paper . In this paper we use only the following result which 
we formulate as 

Theorem 2.3: let H 1 be a finite dimensional Hilbert space (dim H. 1 = n), 
andTL A2 = Tt 1 ATt 1 denotes the 2-particle antisymmetric space generated by 
TL 1 (the 2-fold Grassmann product ofTL 1 ). Let P 2 denote the 1-dimensional 
projection operator onto a 2-particle antisymmetric function g 2 6 TL h2 of 1- 
rank r = 2s possessing the canonical decomposition g 2 = X^=i d I 2i — 1, 2i) 

with 2X1 I & l 2 = 1> where 2*> = >/2^_i A 4>li = ^<fet(02i-i. 4>li) 

is the 2-particle normalized Slater detrminant. Let the identity operator 
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I 1 on Ti 1 possess the decomposition I 1 = Yll=i S + Y17=r+iPi> where 
Pi = l l = <fii ® = 1) ■■■,n) are 1-dimensional mutually orthogo- 
nal projection operators onto the functions \ i) = <pj. Then, the 3-particle 
operator 3P 2 A I 1 possesses the following spectral decomposition 

s=r/2 n 

^ A Jl = E c 1 - 1 & i 2 )( p L-! + + E ^ + • ^( 3P . 2 A Jl )- 

fc=l Z=r+1 

.Here, Pg 2kl , -Pg 2fc (A; = 1, . . . , s = r/2), PjJ (Z = r + 1, . . . , n) are 1-dim pro- 
jectors onto the following functions: 



3 

92k-\ 



7= = , 2 E e<|2»-l,2»,2fe-l) = J - , 3 , g 2 A|2fc-l), 

vi- I Cfc r , V K fe 1 



i = 1 



.3 1 



= -^==== E &|2t-l,2i,2fc) = ^ — __ g 2 A | 2fe)i 

s 

5 f = ^&|2z-l,2*,Z> = V% 2 A|Z), 
i=l 

Ker(3P 2 A I 1 ) denotes the projection operator onto the null- space of the 
operator 3P 2 A J 1 , which is of dimension (g) — n. TTte symbols of the type 
| 2i — l,2i,2k — 1) denote the appropriate 3-particle Slater determinants. 

Now, we can apply results of Theorem 2.3 to find new necessary condi- 
tions for 3-representability according to the formulae of Theorems 2.1 and 
2.2 . The required A max (g 2 ) is either 1- | £ \ 2 min , if the 1-rank of g 2 is 
equal to the dimension n of the 1-particle Hilbert space H 1 (r = 2s = n), 
or A max (g 2 ) = 1, if r is less than n (r < n). Hence, we have 

Theorem 2.4: If D 2 (TrD 2 = 1) is a 2-fermion density matrix, then for 
3-representability it must satisfy for any g 2 G Tt A2 , with the canonical de- 
composition g 2 = Yli=i & I 2« — 1, 2i), Z/te following inequality: 

Tr{D 2 P 2 ) < 1(1- | £ |LJ, (2.3) 
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if the 1-rank of g 2 equals n = dimTL 1 (r = n), 

Tr(D 2 P 2 ) < 1 (2.4) 

if the 1-rank r of g 2 is less than dimTL 1 {r < n). 

In particular, if we take as g 2 the eigenfunctions g 2 of D 2 , we get 

Theorem 2.5: If D 2 has the spectral decomposition D 2 = Yli=i ^^g % w ^h 

9i = EjLi&j I 2j - 1,2.7'), then A; < ±(1- | \ 2 min ), if the eigenfunction 

3 

gf has 1-rank r\ = 2si = n, and Aj < |, if the eigenfunction belonging to Aj 
has 1-rank Vi < n = dimTL 1 . 

It is worthwhile to observe that the eigenfunctions belonging to the max- 
imal eigenvalue A = | cannot have full 1-rank ( r = n = dimTL 1 ) if D 2 is 
3-representable, and on the other hand, the eigenvalues corresponding to 
the eigenfunctions with full 1-rank must be strictly less than |. 

3. Strength of the dual P-condition for fermion 3- 
representability 

It is important to compare the new condition with the already known con- 
ditions for 3-representability because the effort to find A max (g 2 ) for at least 
some g 2 in the general case (arbitrary N) might not pay off. Since the new 
condition is an estimation from the above on the expectation value of D 2 in 
any state g 2 £ 7i A2 , it can be compared with the B- and C- conditions for 
iV-representability of D 2 [2, 3, 21]: 

Tr[B 2 N (g 2 )D 2 } > 0, Tr[C 2 N (g 2 )D 2 } > 0, Vg 2 G Ti A2 , (3.1) 

where 

B 2 N (g 2 ) = I A2 - (AT - 2)L\P 2 A I 1 — (N - l)P 2 G V 2 N , 

C 2 N {g 2 ) = {n-N + 2)L\P 2 A I 1 - (N - l)P 2 G P 2 N . 
In the case N = 3, conditions (3.1) take the form: 

Tr(D 2 P 2 )<^[l-Tr(L\D 2 L\P 2 )l Vg 2 G H* 2 (3.2) 
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Tr(D 2 P 2 ) < '-^Tr{L\D 2 L\P 2 ), Vg 2 e Ti A2 . (3.3) 

In order to compare the inequalities (2.2), (3.2), (3.3), we choose as g 2 
the so called "extreme geminal" [2] g 2 xtr of 1-rank r = 2s = n = dimH 1 
possessing the following canonical decomposition: 

n/2 

9extr = ^2]/-\ 2i-l,2»). 

The required by inequality (2.3) | £ \ 2 n in = n' an< ^ therefore the new condition 
gives 

Tr(D 2 PlJ < 1(1 - I). (3.4) 

On the other hand, since L^P 2 ^ = £? =1 \P} > P i =1 *><« I. £?=i ^ = ^, 
we have from both (3.2) and (3.3) the same result 

Tr(D 2 P 2 ex J < \{l - i). (3.5) 

Comparing (3.4) with (3.5), we see that the new condition is stronger than 
the B and C ones. So far, we were unable to establish the relation between 
the new condition (2.3) and the G-condition [13] for 3-represent ability. This 
would be important because it is known [11, 20] that the G-condition implies 
the B- and C- ones. 



4. The strengthened B-condition 

While comparing the new condition with the B-condition we found out that 
the B-condition could in principle be improved for N odd. B 2 ^(g 2 ) £ Pjj, 
Vg 2 e H A2 means that the operator (^) Bff (g 2 ) A / A ( JV_2 ) > 0, Vg 2 G H A2 . If 
^mfn * s t ne minimal eigenvalue of the operator (^Bfj^g 2 ) A / A ( Ar_2 ), then 
the operator ( 1 ^)B 2 v (g 2 ) A I A ( N ~ 2 ) - A^£l AN is positive-semidefinite and 
therefore 

{^)B 2 N {g 2 )-K B ±I A2 eV 2 N , (4-1) 

giving a necessary condition for iV-representability (the strengthened B- 
condition). It is known [3] that for N even = 0, Mg 2 £ H A2 . For 
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N = 3 and g 2 xtr , the eigenvalue A.^^ xtr ^ is greater than zero, and can be 
found explicitly. The operator 

W 2 (g 2 extr ) A I 1 = 3(1 - -V 3 - 2(3Pl tr A I 1 ) 

possesses the following minimal eigenvalue 

A2r r) =3(l-i)-2(l-^) = l + i 

where Theorem 2.3 has been used. Hence, the strengthened B-condition in 
this case is 

ZBM xtr ) - aST^ 2 = ^^V 2 - 2(3i*„ J G f* N , (4-2) 
and therefore 

Tr(D*PlJ < \(1 - \). (4.3) 

Comparing inequalities (3.4) and (4.3) we see that both the new conditions 
for 3-representability in the case under consideration (g 2 = g 2 xtr ) give the 
same bound from above on the expectation value Tr(D 2 P 2 extr ). 

It seems that the above results suggest, it would be worthwhile to make 
an effort to extend the new conditions to arbitrary N for at least some g 2 for 
which the maximum eigenvalue A max (g 2 ) of the operator (J%)Pg A I A ( N ~ 2 ) 
can be found, e.g. the extreme geminal g 2 xtr - So far, we have obtained 
only partial information about the spectral decomposition of the operator 
i^Pgextr ^ I A ( N ~ 2 ) for arbitrary N. However, at least in this case it is 
realistic to succeed. The other known conditions for iV-representability could 
be treated in a similar way as the P- and B- conditions considered in this 
paper, provided the maximal and minimal eigenvalues of the appropriate 
operators were known. The obtained results concerning arbitrary N will be 
published later. 
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